Shell sources as a probe of relativistic effects in neutron star models 
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A perturbing shell is introduced as a device for studying the excitation of fluid motions in 
relativistic stellar models. We show that this approach allows a reasonably clean separation of 
radiation from the shell and from fluid motions in the star, and provides broad flexibility in the 
location and timescale of perturbations driving the fluid motions. With this model we compare the 
relativistic and Newtonian results for the generation of even parity gravitational waves from constant 
density models. Our results suggest that relativistic effects will not be important in computations 
of the gravitational emission except possibly in the case of excitation of the neutron star on very 
short time scales. 
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I. INTRODUCTION AND OVERVIEW 



For most astrophysical objects Newton's classical theory of gravity gives a fully satisfactory description. Only 
when gravitational fields become strong need one consider the possibility that general relativistic effects may play a 
significant role. A standard index of field strength is GAI/Rc 2 , where M is an object's mass, and R is its characteristic 
(N size. This index is of order unity for black holes and for the universe itself, and much smaller than unity for almost all 
stars, galaxies, and other astronomical entities. One exception is neutron stars; for compact neutron stars GM/Rc 2 is 
on the order of 0.2, and exotic equations of state could lead to even larger values. Despite this, Newtonian gravity is 
■ used almost exclusively in studying neutron stars. The obvious reason is the significant increase in difficulty in giving 
I a fully relativistic treatment of neutron star structure, and the enormous increase of difficulty in dealing with fullly 
. relativistic dynamics, i.e. , with oscillations of neutron stars. 

Newtonian physics has been used even in studies (see for example Refs. of neutron stars as sources of 

gravitational waves. Of course, Newtonian gravity per se has no gravitational waves, so the typical procedure is 
to compute gravitational wave generation as a postprocessing step. More specifically, Newtonian gravity is used to 
find the fluid motions inside a neutron star associated with some event (core collapse in a supernova, precession of 
O a rotating neutron star, etc.). Those fluid motions are then used as sources in the "quadrupole formula" of general 
relativity, just as known charge motions would be used in the dipole formula of electromagnetism. If Newtonian theory 
5h . predicts periodic oscillations, the "adiabatic approximation" can be used to give the damping of the oscillations due 
^P" to gravitational wave emission: the energy of the oscillations is taken to decrease at the rate at which gravitational 
^ , waves remove energy. 

The sufficiency of this approximate procedure caused little worry until relativistic modes of oscillation of neutron 
stars were discovered Ji]^] that had no counterpart in Newtonian theory. The so-called w modes are qualitatively 
oscillations of the spacetime, like black hole quasinormal (QN) modes, rather than oscillations of the neutron star 
material like the / (fluid) and p (pressure) modes of the Newtonian description. Partly due to the existence of these 
w modes, the general question of the sufficiency of Newtonian theory for dealing with gravitational wave processes 
was emphasized by Andersson and Kokkotas 0, and stirred considerable interest. 

To give a clear answer to this question requires a specific and astrophysically plausible event for which gravitational 
wave generation can be computed in both Newtonian theory and fully relativistically. One would want, for instance, 
initial data for the fluid and spacetime of a neutron star formed in a supernova core collapse, but the possibililty 
of giving such initial data is at least several years off. A more tractable model was needed. Three separate groups 
IpHlOl studied the problem of emission of gravitational waves by a relativistic neutron star due to the close passage 
of a perturbing particle. This model had the advantage of definitiveness; there was no freedom in choosing (and 
biasing) the initial spacetime perturbations to be particularly larger or smaller than the fluid perturbations. This 
model allowed an investigation of whether the excitation of w modes was significant, but the model had three serious 
shortcomings, (i) It was too restrictive. The only parameters were the two constants specifying the particle orbit (say, 
energy at infinity and impact parameter). There could be significant excitation only if the particle passed close to the 
neutron star, and this constrained the perturbation to be neither very fast nor very slow. Not only was the timescale 
limited, but it was coupled to the choice of location of the perturbation, (ii) It was difficult to separate the radiation 
due to the neutron star (the radiation of interest) from the radiation coming (in some sense) from the orbiting particle, 
(iii) It was not clear how to compare the fully relativistic computation from the Newtonian computation; at least 
there was no attempt to do this. 
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We introduce here a very different model for investigating the importance of relativistic effects in neutron stars, 
and possibly for answering other questions. We consider a perturbative spherical shell around a neutron star at some 
radius i? s hcii- In that shell we dictate the time dependence of a multipole of surface mass-energy density of the shell. 
There is no equation of state of the shell material constraining our choice. The equations of motion of the shell 
fix the surface stress once we have specified the surface energy density, so picking that single function of time fully 
specifies the source. In this manner we can independently choose where the perturbation of the neutron star arises 
and what its timescale is. We can probe the response of the star to close and far perturbations with slow motions 
or fast motions. The shell probe has the nice feature that it is straightforward to do the calculation in Newtonian 
theory, so that a comparison can be made with the relativistic result. A further advantage is that both the Newtonian 
and relativistic computations allow a separation of the gravitational waves from the shell and those from the star. In 
the Newtonian computation, this is completely straightforward. The waves from the star are those found from the 
quadrupole formula applied to the motions only of the stellar fluid. In the relativistic calculation the separation is 
only approximate, and more care is needed. From the star+shell results it is necessary to subtract the radiation from 
the shell itself generated in the background spacetime of the star. (The details of this procedure will be given below.) 

Our main purpose in the present paper is to present the method of using a shell probe and to display its advantages. 
For that reason we limit the application of this method to the simplest model of a star, the homogeneous incompressible 
perfect fluid (HIF) model. (See 11 for details about the QN modes of this model.) For this model the only mode 
associated with motions of the stellar material is a single / mode. An important element of our results is that we 
will present answers not only about excitation of w modes, but about the difference in the Newtonian and relativistic 
predictions of excitation of the / mode. In the interest of brevity we limit the analysis to even parity perturbations. 
The excitation of w modes for odd parity should not be remarkably different from the excitation in even parity, and 
odd parity motions do not couple to fluid motions. 

The paper is organized as follows: The shell source model is introduced in Sec.||, and the equations governing even 
parity perturbations due to the shell source are given in Sec. HI. Some details of the computational implementation are 
given in Sec. IV along with a discussion of the method used for subtracting the shell contribution from the relativistic 
calculation of the gravitational wave due to the shell and star. Numerical results are presented and discussed in Sec. 0, 
and conclusions are given in Sec.^. Details of the Newtonian calculation are given in the Appendix. Throughout 
the paper we use geometric units G = c = 1, the metric signature ( — h H — h) and the conventions of Misner, Thorne 
and Wheeler @. 



II. MODEL: PERTURBATION OF A STATIC STAR BY MATTER MOVING ON A SPHERICAL SHELL 



We start with a static and spherically symmetric spacetime background metric 

ds 2 = -~e u{r) dt 2 + e x ^dr 2 + r 2 [d9 2 + sin 2 9d<p 2 } , (1) 

describing both the interior and exterior of a star of a barotropic, ideal fluid of mass M and radius r = R. The stress 
energy of the fluid is 

T a/ 3 = (p + p)u a up + pg af3 , (2) 

where p is the mass-energy density and p = p(p) is the pressure. The mass function m(r) is defined by e~ x ^ — 
1 — 2m(r) /r, and the structure of the stellar interior is found by solving the hydrostatic equilibrium equations of 
general relativity. (See, e.g. , Eq. (3) of |lC|].) For simplicity, we limit considerations to homogeneous incompressible 
fluid (HIF) stellar models whose unperturbed interior metric is given by Eqs. (5) and (6) of jlCfl . The exterior metric 
is simply the Schwarzschild metric, with m(r) equal to a constant M, and v(r) — —X(r). A spherical thin shell of 
coordinate radius r = i? s hcii > R surrounds the star. We treat the shell as a perturbation of the spacetime inside and 
outside the star and we analyze the perturbations only to first order in the parameter of the perturbation. In this 
order, the shell has spherical geometry described by the 3-metric, 

ds\ hell = - ( 1 - dt 2 + R 2 hcll [d9 2 + sin 2 0dip 2 ] , (3) 

\ -Kshcll / 

induced by the Schwarzschild metric. The metric of the perturbed spacetime can be written as 

9a/3 = 9„p + h »/3 (4) 

where the "(0)" index denotes the background solution, that of Eq. ([!]). The Einstein equations to first order in 
perturbations are 
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SG a ^ = 8n[6Tj Rnid + ST a i en }. (5) 

The perturbed stress energy has two contributions. One, denoted by 5T a g uid , is that of the fluid star perturbations 
and is nonzero only inside the star. The other, denoted by 6T a f h u , is the stress energy of the matter in the thin shell 
and is nonzero only outside the star. Its form, in the coordinates of Eq. (^), is 

ST^ ell = y/1 - 2M/r S a ?6(r - R shcll ) , (6) 



where 



S a ? = lim / rpajj ar _ 

is the surface stress energy of the shell. (See e.g, 

A. Even parity equations of motion of the matter in the shell 

From Eq. (|J), we obtain the equations of motion of the shell, ST af3 s hcii:/3 = 0, where ; denotes the covariant 
derivative in the Schwarzschild spacetime. Taking Eq. (^|) into account leads to the restriction S ra — and to the 
partial differential equations, 

S ab lb = 0, S ra - a = 0, a,b = t,6,tp (8) 

for the components of the surface stress energy tensor. Here "i" is the covariant derivative with respect to the shell 
3-metric in Eq. (||) . 

Due to the spherical symmetry of the shell, we can decompose Sqq, Soi, Sij]i, j = 6, (p in scalar, vector and tensor 
spherical harmonics respectively. Restricting attention to the even parity harmonics, we write these decompositions 
as: 

S 00 = J2 S iomo(0) (9a) 
Soe = J2 S to(t)§QYio(9) (9b) 

Sa = !>«(*)*«>« (0) + Sf (t)* mj (6)] . (9c) 
i 

Here Y w is the scalar spherical harmonic and $jo w / sin 2 9 = § me = Yi ,^>i ogv> = and ^wse = d 2 /d9 2 Y m , ^i 0vv = 
sm 9 cos 6d / dOYiQ, ^loetp = are the even parity Regge- Wheeler |l3j tensor harmonics. Azimuthal symmetry guaran- 
tees that the azimuthal index m does not enter into any of the equations after multipole decomposition, so with no 
loss of generality we consider only m = axially symmetric motions. 
Upon substitution of Eq. (^) into Eq. (||) we get, 

Rshen ^^ + 1)^ = (10a) 



1 - 


2M/iUell 


dt 








1 - 


2M/R sheQ 


dt 






\9 D 10 



?fo + W + 1) - 1] Sf = (10b) 

2Sf + 1(1 + l)S? = . (10c) 

As these equations show, we only have one degree of freedom. The choice of the surface mass-energy density Sf (t) 
uniquely determines all the other components of the shell's stress energy through Eq. ([To|). In this work, we make the 
choice 

-at 2 

Sf (t) = e— , (11) 

where e is the perturbation parameter. Since all perturbation equations will be proportional to e we will omit it 
henceforth. The use of a Gaussian time dependence for the surface density on the shell gives us a source that it 
localized in time and allows us a choice of timescale for the process that drives the stellar fluid motions. 
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III. EQUATIONS GOVERNING EVEN PARITY PERTURBATIONS 



A multipole decomposition of the even parity quantities in Eq. (|5|) leads to a set of coupled partial differential 
equations in the variables t, r for the coefficients of the metric perturbation h a p and velocity of the fluid star 
p5[ . We adopt the notation of Regge and Wheeler pi, Thorne and Campolattaro Q and of Moncrief |l6[ . For 
simplicity, we make the Regge- Wheeler JT^] gauge choice, which for even parity means that the only nonvanishing 
metric perturbations, for a particular are 

h m = e^H l \r,t)Y m (6) (12) 

h Or = H l 1 o (r,t)Y lo (0) (13) 

h rr = e x ^H l 2 °(r,t)Y l0 (9) (14) 

h ik =r 2 K m {r,t)^ l03k , j,k = 9,<p, (15) 

where $;ojfc is one of the even parity tensor harmonics defined above, after Eq. ([)]). It is useful to divide the perturbed 
Einstein equations (||) into those governing perturbations inside the star and those for perturbations outside the star. 

A. The interior equations 

Inside the star the only nonzero stress energy is due to the perturbed fluid. Its independent components, are 

ST ° uid = -dp (16a) 

ST r fluid = ST e 9 Ruid = 6T V £ uid = Sp (16b) 

5T g uid = (p + p)u Su a , (16c) 

ST a fluid = (p+p)Su a u°, a = r,9,cp, (16d) 

where Sp : Sp are the Eulerian changes in density and pressure, u° = e~ vt ^ r ^ 2 is the only nonzero component of the 
velocity of the unperturbed star and Su a is the velocity of the perturbed fluid. It is convenient to introduce, at this 
point, the quantity 

thmJlL.. (17) 
p + P 

For barotropic fluids, Sh is the Eulerian perturbation of the relativistic enthalpy and 

Sp = { P±ffsh , (18) 
PI 



where 7 is the adiabatic index, 



^E±£fl±. (19) 
p dp/ar 



We decompose the stress energy components in Eqs. (|l^) into spherical harmonics, and for a single multipole have 

5h = 5hi(r,t)Y l0 (9) (20) 
5u° = ^e-'^/ 2 H l o °(r,t)Y la (0) (21) 

-[u(r)+\(r)]/2 a 

Su r = g Wlo (r,t)Y m (0) (22) 

r z at 



6u 6 = ^^Q-V l0 (r,t)-Y lQ (9). ,2.-!) 



With a similar decomposition of the perturbed Einstein tensor into tensor harmonics, Eq. (^) leads to a set of coupled 
equations for Hq°, H[°, K m , Hlf, Wiq, Viq and ShiQ. One of the equations |T^] provides the important simplification, 



H l °(r,t)=H l °(r,t) . (24) 
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Using this result, the other equations, can be reduced to a coupled system of two equations |17fl in which all fluid 
functions have been eliminated and the only dependent variables are the metric functions Hq°,K 10 . For barotropic 
fluids, these equations are 

e [A<r)-Mr)]/2[ e [Mr)-A<r)]/2j^ ir + 2 (-1 + 1^ Kf - e X ^-^K% -(I- l)(l + 2)^K W - H l ° rr 
+ (2 + \ L _^r_\ H lO PJJ+}l e Hr) _ 2 (1 _ rVr) e[ A(r)-3Kr)]/2 [e [3Kr)-A(r)]/ 3l/ ] ^ 

\r 2 2 J ' [ r z r z 
-8ne x ^ (p + p) } i?Q + e^-"^H l % = , (25) 



and 



K 



/(i 



r 2 j»7 

P + P \ ri 0,r 



/o 



■(i - + 2) f 1 + £ + A''° - e A(r)-(r)£±£ lr 

V P7 / P7 



P7 / r 



2r 2 



pry \ r 



2r 2 



r j»7 



M°=0 



(26) 



For HIF models, the adiabatic index 7 is effectively infinite since p =constant, and thus all the terms proportional to 
(p + p)/p r y in Eq. (p6|) are zero. 



B. Reduction of the exterior equations to the Zerilli equation 



In the Schwarzschild spacetime outside the star, it can be shown |15| |18| that all the perturbation equations can 
be obtained from the two first order equations 



1(1 + 1)H[° + 2rH l ° t - 2r 2 K% + ? M n J* K* = D l0 (r, t) 



1 - 2M/r 



2M 



=^ H {° +(!- — ) HZ r K" = B l0 (r, t) 



2M\ 



10 TjlO 



r I 



together with the algebraic identity 



F 



(/-!)(/ + 2) + 



6M 



2r 2 



K l ° + 2rH 
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ttIO 
H 0,t 



M 



(/-!)(/ + 2) + 



2M(r ~ 3M ) 



r(r - 2M) 



K 



10 



1(1 + 1) 



Hi - 



l-2M/r"' m ^'" x Mt 
and the shell's equations of motion, Eq. (|io|). The source terms in Eqs. (|27|), ( |2S| ) are 



D l0 = 32nr^Jl-—6(r-R shcn ) 
at V r 



B w = 16tt { S] 



2M 

1 S (r - Rshen) 

r 



Equations @, (|§ and 



for the Zerilli Ifl8fl function, 



/i; dt 

in turn, can be combined into the single wave equation, 
d 2 Z l0 d 2 Z l0 



dt 2 



,.;,,:•> ■ V,(r)Z w =5 w (r,t) 



Zm(r,t) 



r(r - 2M) 



/Hi 



(rA + 3M)(A + l) 
In Eq. (B2), r* is the usual tortoise coordinate, 



[H£ - rK lu r 



A + l 



-A" 



(27) 
(28) 



(29) 

(30) 
(31) 

(32) 
(33) 
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the constant A is 



and the potential Vi is 



r* = r + 2Mlog[r/2M - 1] + constant . 

A _ (i-m + 2) 



2A 2 (A + l)r 3 + 6 A 2 Mr 2 + 18AM 2 r + 18M 3 / 2M 
' W ~ r 3 (Ar + 3M) 2 



The source term in Eq. is 



ro( ' } = — W) — ~ 



C i0 (r, t)*[r- J2 



shell 



h(r) 



16tt / 2M 



ft 



shell 



ft 



hell 



3/2 



Sf (t)*(r- ftshell) 



where 



C l0 (r,t) 

m 

h{r) 



16nRl 



hell 



l(l + l) y /l-2M/R Bbfia 
6M 2 + 3AMr + r 2 A(A + 1) 

r 2 (3M + Ar) 
3M + Ar 
r-2M ' 



S? (t) 



(34) 



(35) 



(36) 



(37) 

(38) 

(39) 
(40) 



Together with appropriate boundary conditions, the set of equations fl25[), ( |2q ) and (|32|) govern the even parity 
perturbations of the star due to the shell. 



C. Boundary conditions 

1. Regularity at the center of the star 



The system of equations (|25-26) admits two linearly independent regular solutions. Near the center they admit the 
series expansion [ p7| , 



K l0 {r, t) = k (t)r l + k 2 (t)r l + 2 + 0{r l+i ) 
H l °(r, t) = k (t)r l + h 2 (t)r l + 2 + 0{r l+i ) 



(41a) 
(41b) 



The r-independent "constants," ko, k 2 , h 2 are a priori arbitrary. When Eq. ( pq) is expanded about r = and the 
above expansions are used, we find 



rs i 2 ' 




\l 2 11, 1 




h 2 + 


2 + T* + 9 









k 2 - y po [21 + l 2 -3] k Q = 



(42) 



where po,po, 7o are the values of the density and pressure at the center of the star and we used the fact that 7 = 00 
for a HIF. Equation ( f42|) allows us to eliminate fco, so that the only remaining unknown coefficients are k 2 and h 2 . 
These turn out to be fixed, up to an overall scaling, by conditions at the stellar surface. 



2. Vanishing of the pressure at the stellar surface 

The pressure must vanish at the perturbed surface of the star by definition of that boundary. This amounts to the 
vanishing of the Lagrangian perturbation of the pressure at the stellar surface r = R, 



G 



Sp(R) 



-\(R)/2 



R 2 



-W(R,t)p, r (R) = 



(43) 



When the density p vanishes at r = R, the second term in Eq. (|4^) is zero and this condition is the same as the 
vanishing of the Eulerian perturbation of the pressure. But this is not so when p(R) ^ as is the case for a HIF. 

After differentiating Eq. (43) twice with respect to time and doing a decomposition in spherical harmonics, we can 
rewrite it purely in terms of the metric functions K l0 ,Hjf and its derivatives [fl7| . The resulting expression, when 
p(R) / > 



-K 



ID 



-H, 



O.rtt 



V Mr) K m 1 (l-m + 2) pA(r) V u , 

4 r 2 V 2 r 2 



,. v , K iu _ \{r)-u(r) K h 
e tt e 



WO 



1 K 

27 



4r 2 



2r 



2r 2 



H 0,tt + KKr) e ' —, 2 H 



',_+_ 
4r 2 



11 .tt 



r=R 



(44) 



3. The relation between the interior and the exterior metric functions 

If the density of the star is not zero at r = R, the first radial derivative of K 10 ,Hq° is not continuous at r = 
R, although the functions themselves are. Denoting by a superscript "+" the exterior metric functions and by a 
superscript "— " the interior ones, we have pM, for a particular multipole I, that 



K+(R,t) = K~(R,t) 
H+(R,t)=Ho(R,t) 



K+(R,t) = K~ tt (R,t) 



2M 



2R 2 



1(1 + 1) [1 ff - - 2rH 



r J 

1(1 + 1) 
2M y ' Hq 



2r 



r-3M 



r-2M 



■K 



,rU 



r=R 



(45a) 
(45b) 



(45c) 



If p(R) = then both metric functions and their first radial derivatives would be continuous at r = R. 



D. Fourier transform of the perturbation equations 



The simplest way to solve the above partial differential equations in r, t is to write all time dependent quantities as 
Fourier integrals, reducing the problem to that of ordinary differential equations in r. Thus we write Zio(r,t) as, 

1 f°° 

Zlo(r,t) = — e- lut Z w (r,io)dw (46) 



2?r . 

transforming Eq. (|3^ ) into a second order equation for Zio(r,u>) 

d 2 Z l0 



dr* 2 

where 



+ [uj 2 -Vi(r)]Z l0 = -S w {r,Lu) (47) 



S l0 (r,io) = / S w (r,t)e iuJt dt (48) 



is the Fourier transform of the source term. We can proceed similarly with Eqs. (25) and (|26|). The resulting equations 
for -ffo°( r i w )> K l0 (r, lu) can be obtained directly from Eqs. ( p5| ) and ( |26| ) by substituting — iu> for dt, and by replacing 
functions of r, t by their Fourier transforms. 

At infinity, we impose the boundary condition that the wave be purely outgoing 

Zi (r^aa,u=t-r*) = ±- [ A w (cu)e~^ u duj . (49) 
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The even parity gravitational energy, radiated in a single multipole component, is then given |l9[ by 



dui 



1 (l + 2)\ 
64tt 2 (l-2)\ 



w 2 \A l0 (u)\ s 



(50) 



From the wave function Ziq(t, t) we can construct the multipoles of the metric perturbations. The transverse traceless 
(TT) part of the metric perturbation is of particular interest, since it completely characterizes the radiation at infinity. 
Extracting the TT part is easiest if we realize that the Zerilli function Zio(r,t) is numerically equal to Moncrief's |lq | 
gauge invariant wavefunction. Since the Moncrief invariant can be evaluated in any gauge, we choose the gauge to be 
asymptotically flat so that all multipole perturbations fall off faster than 1/r except the TT components. From this 
procedure we find that the TT perturbations are related to the Zerilli function by 



1 = 2 



2(1-2)1 jk 



(51) 



where 



10 jk 



1(1 + 1) 



4> 



10 jk 



(52) 



is the orthonormal even parity TT tensor harmonic and ^iq, $zo are the tensor harmonics introduced in Eq. 



IV. COMPUTATIONAL IMPLEMENTATION 



A. Solution for Aio(lo) 

If the background spacetime is due to a star, a solution of Eq . ([47| ) must be found that corresponds to outgoing waves 
at infinity and that matches the regular solution of Eq. (|2q), (|26|) at the unperturbed surface of the star according to 
the junction conditions of Eq. (|45"|). The Green function solution is found in the usual way. (See e.g. , Ref. (l(J.) We 
define y° ut (r,uo) as the homogeneous solution of Eq. ( |47| ) with the asymptotic form 

!/?"*(»•>«) -» e lur " r*->oo. (53) 



For our second independent solution of Eq. ( |47| ) we start by finding, in the stellar interior, a solution of Eqs. fl25j), ( |26| ) 
satisfying the condition Eq. ( pl| ) at the stellar center. Our second solution y\ e9 (r,u>) is taken to be the homogeneous 
solution of Eq. (|47|) that joins to the interior solution through the matching conditons of Eq. (ff5j). 
We then define the Wronskian of these two homogeneous solutions, an r independent quantity, to be 

rfyout j reg 

^)-vV 9 %r-yr%r- (54) 

With the above definitions, and from the Green function solution, we obtain (see, e.g. ,Ref. jl(|) the Fourier amplitude 
Aio defined in Eq. @ 

i /<oo reg / \ 

A l0 (u) = -—— / S l0 (r,u) Vl I; dr . (55) 

Combining Eq. (jnj, Q, @ and @ we get explicitly, 



yl e9 (R she ii,ui) 
l-2M/R shcU 



(Z — 1)(Z + 2) V Ml(l + l)J Mh(r)l(l + 1) 



(56) 



What is required for a solution, then, is a numerical determination of y i res (i? s i lc ii, w) and its radial derivative. 
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B. Numerical method to find y rea for a HIF 



The numerical pro blem of finding j/; re9 (i? s hoii, and its derivative for a HIF, can be divided in two parts: Integration 
of Eqs. (p5| ) and (|2^) from the center of the star to r — R and integration of Eq. ( [t7|) from r = R to r = i? s hcii, so 
that we can evaluate y\ e9 , dy/ eg /dr at this radius. 

To find K~(R, u), Hq(R, lu) and K~(R, ui), for a HIF we first find an "A solution" by starting from the center with 
with fco = 1, &2 = and /12 obtained from Eq. (|42|); we then integrate Eqs. (|25|), ( p6| ) out to r = R, to find the "A 
solution" there. The procedure for the "B solution" at r = R is the same except we start with the central conditions 
fco = 0, fc2 = 1. The general solution for and K~ can be written 

H f J = a H§~ + (3 R- = a K A ~ + K B ~ . (57) 

where a and /3 are arbitrary constants. The overall scale of y ; re9 (r, lj) is arbitrary. Since yj eg (r,u>) occurs both in 
the integrand in Eq. ( |55| ) and in the Wronskian in the denominator, the scale cancels out so only the ratio a/(i is of 
importance. This ratio can be found by substituting the expression in Eq. ( |57| ) and its derivatives in the vanishing of 
the Lagrangian pressure condition Eq. (]44|) . 

We can next compute K + (R,u), H + (R,ui) and K+(R,ui) using Eq. (^). The final step is to use Eq. (|3|) to find 
the starting values for integrating the Zerilli equation, 



dr 



3M 



( J RA + 3M)(A + 1) 
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and to integrate the Zerilli equation Eq. (|32| ) out to r = i? s hcii in order to find 2/; re9 (i? s hoii, w) and dy^ 3 / dr{R s \ lc \\ 



C. Radiation due only to the shell 

The waveform Z[q, its Fourier transform Aio(lo), and the energy computed in Eq. (|5^) refer, of course, to radiation 
from the star and the shell. As might be expected the radiation that can be attributed directly to the shell is 
much larger than the radiation from the perturbed fluid motions in the neutron star. In order to have the clearest 
comparison of Newtonian and relativistic predictions of radiation from the neutron star, it is useful to remove the 
contribution due to the shell. 

For the Newtonian computation presented in the appendix, this presents no problem. As is evide nt in Eq. ( Al§| ), 



the density perturbations due to the shell and due to the star are distinct, and it is evident in Eq. ( |A2l| ) that their 
contribution to the quadrupole moment are distinct. This clear distinction does not exist in the relativistic calculation. 
The equations of Sec. Ill contain information about the shell entangled with information about the star. To get an 
approximate idea of what radiation can be ascribed to the star, and not to the shell, one can compute the waveform 
due (in a sense) to the shell itself, and subtract this waveform from the total star+shell waveform. This subtraction, 
however, is somewhat subtle. 

In particular, it is not useful to consider the shell in a flat spacetime background. If we consider a shell with the 
surface stress energy of Eqs. @, (|ll|) radiating in a flat background, and the same (in some sense) shell, and surface 
stresses, radiating in another fixed background, there will be a large difference simply due to the different radial null 
geodesies, the spacetime lines on which perturbations propagate. If we are to have the same shell-only radiation 
as that due to the shell in the shcll+star problem, we must have the two signals propagate on the same spacetime 
background. To accomplish this, we describe the waves at radius r > R for the shell-only problem with the same 
Zerilli equation as we use in the shell+star problem. Equation (B2|) is then part of both problems for r > R. In the 



shell+star problem the perturbations in the interior are, of course, treated with the equations of Sec. [II. For the 
shell-only problem, we need instead something like a Zerilli equation suitable to the fixed spacetime of the stellar 
interior. Such an equation is provided by the equation for the propagation of massless scalar perturbations in the 
neutron star background. These waves do not excite oscillations either of the spacetime or of the neutron star fluid. 
They are described only by an equation identical in form to Eq. (|3^), but with no source, and with the Zerilli potential 
of Eq. (|36|) replaced by 
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scalar 
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(60) 



(The idea of freezing the fluid perturbations to study the w modes described in |2(J is somewhat similar in spirit to 
our method but is different in practice.) With this mixed mathematical description we find the regular solution of the 
Zerilli equation and then compute, using Eq. (56), the Fourier transform A;q c11 and, with Eq. ([49|), a wave function 
•ZJq 1c11 . By subtracting these, respectively from the Aio and Ziq for the shell+star, we arrive at results Ajg ar , Zf^ 1 
meant to describe radiation only due to the oscillations of the neutron star fluid. From the square of ^4jo ar , we can 
compute an energy spectrum and a total radiation energy attributed to the motion of the stellar fluid. 

Our method of finding the radiation from the shell alone is, of course, only an approximation. Einstein's equations 
couple oscillations of the fluid and oscillations of the spacetime, so that there can be no completely meaningful way 
of separating the two. One sign of this is that the potential Eq. ( |60|) is not unique. The potential does not influence 
the radial null geodesies, and any potential with the same general behavior at r — > 0, is equally good. In situations 
for which the details of the potential are important, our method of computing shell-only radiation is not justified. 
But such details are not relevant for must of our models. This can be seen in the reasonable success of our method 
in removing appearances of shell radiation in the star-only results presented below. 

An additional important point to understand about our shell subtraction method is that in principle it should 
subtract the w modes. Since the w modes are due to the spacetime background, not to the fluid motion, the shell-only 
radiation should have w modes, and the star-only waves that result from subtraction should have only waves due to 
fluid motion. This will be discussed further in connection with the numerical results presented below. 



V. NUMERICAL RESULTS 



We focus first on a very compact HIF stellar model of radius R = 2.5M. By starting with an extreme, though 
astrophysically implausible, model we will be able to see relativistic features that will be absent in less compact, 
and more plausible, models. In Fig. ^(a) we present the waveform Z20, the Zerilli function for quadrupole (I = 2) 
perturbations in the case that the Gaussian parameter in Eq. (|ll| ) is a = 0.1Af~ 2 . The shell was placed at a large 
distance from the star, i? s hcii = HOM, to have the waveform clearly display the time profile of events. The first burst, 
at around u = t — r* = — 160M is radiation coming directly from the stress energy of the shell. Later, at around 
u = t — r* ~ lOOAf , a burst arrives representing the ingoing radiation from the shell that has "reflected" off r = 0. 
At around the same time, radiation from oscillations of the fluid and central region of the spacetime arrive. 

Two curves are shown in Fig. [j](a). One is the waveform for the r elativistic star+shcll. The second curve is 
that for the shell itself, computed by the method described in Sec. IV C. The two curves are nearly identical up to 
around u = 160Af, confirming that the early radiation is that due to the shell. From u — 160Af to around 400Af or 
more, there are damped oscillations of the dominant (least damped) w mode. The w modes depend on the details of 
wave propagation in the innermost strong field regions, and are not the same for the star+shell problem and for the 
somewhat ad hoc spacetime we constructed for the shell-only problem. The w mod e frequ ency for the constant density 
R = 2.5M model is to = (0.42 + i0.0217)A/-\ while that for the spacetime of Sec. |lV C| , w = (0.438 + i0.029)M~ 1 , is 
more rapidly damped. We would, of course, find a different w mode frequency if in the shell-only problem we used a 
potential other than that in Eq. (p0|). 

It is clear that there is no point in subtracting the shell-only radiation from the star+shell. The difference waveform 
would contain an artifact corresponding to the difference of the w modes, and hence would be dominated by features 
completely irrelevant to fluid motions of the neutron star. Figure [l] (a) helps to demonstrate that subtraction is 
pointless when w mode radiation is of importance in the waveform. It is only the fact that w mode radiation is a 
minor feature for realistically compact stars that makes subtraction useful. 

In addition to the difference in frequency of the w modes there is another, more important, difference between the 
two curves in Figure |l|(a). The star+shcll curve shows an oscillation with imperceptible damping at a frequency less 
than half that of the w mode. This is the / mode of the fluid of the neutron star. This mode is missing, as it should 
be, from the shell-only computation. 

Both curves in Fig. |l|(a) show relativistic results. A comparison between the Newtonian and relativistic star+shcll 
results would not be of much use. The slow motion condition underlying the Newtonian approximation would be 
strongly violated since the radius of the shell i? s hoii = HOM is much larger than the characteristic time scale of the 
shell stress-energy oscillations (At ~ a -1 / 2 ~ several M ) or of the modes of oscillation of the spacetime or of the 
stellar fluid. The Newtonian computation, based on the slow motion approximation (time scale <C light travel time 
across source) would be completely inappropriate. In order to construct a more justifiable Newtonian comparison 
for a R = 2.5M HIF model, we consider in Fig. [l](b), a shell with radius i? s hoii = 5M and with Gaussian parameter 
a = 0.001AT -2 (and hence timescale ~ 30Af). The waveform in this case is not of primary interest, since the radiation 
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from the shell and from the stellar fluid will not be clearly distinguishable. The energy spectrum, however, gives the 
answer to the most important questions, ft shows, for example, that there is negligible difference between the broad 
spectra of the Newtonian and the relativistic results. But the broad spectrum is due to the shell. Of more astrophysical 
interest is the / mode excitation. The relativistic computation shows a lower frequency / mode containing almost an 
order of magnitude more radiation energy than the / mode peak of the Newtonian computation. This conclusion, of 
course, applies to a model that is too compact to be astrophysically relevant. 

The feature in the relativistic spectrum at u> ~ O.SM^ 1 is due to the repetition of the shell radiation with a 
time delay of At ~ 10AT. This produces a modulation of the form cos 2 (wAi/2), and hence a dramatic decrease at 
w ~ it/ At ~ 0.3. The location of such features is dependent on i? s h n and, of course, is unrelated to the physics of the 
neutron star. It is worthwhile noting that these features are absent in the Newtonian spectrum. Since that spectrum 
is based on the slow motion approximation the whole star+shell source is treated as if it radiates in phase, and there 
can be no repetition of the shell radiation. 

In Fig. U we show results for a HIF star of radius R = 5M, a typical radius of a neutron star. The shell is located at 
-Rsheii = 10M and the Gaussian parameter for the time profile of the shell stress energy is a = 0.01A/~ 2 . In Fig. g(a) 
the dotted curve gives the full Zcrilli quadrupole waveform of the star+shell computation, and shows why subtraction 
of shell radiation is useful: the waveform is completely dominated by the shell radiation, and the much smaller / 
mode is barely visible. The solid curve is the waveform from the shell-only computation. It is clear in Fig. ||(a) that 
the two waveforms are nearly identical at early times, and at late times are different in that the star+shell result 
has / mode oscillations, while the shell-only waveform, of course, does not. This suggests that subtraction will be 
very effective in isolating the radiation due to the stellar fluid, and this turns out to be true. The star-only curve 
(the result of subtracting the shell-only from the star+shell) is given in the inset to Fig. ||(a) and shows the / mode 
oscillation as a dominant feature. 

Figure ||(b) shows the energy spectrum for the subtracted (i.e. , star-only) case shown in the inset of Fig. ||(a) . Both 
the relativistic and Newtonian spectra are given. The two spectra are roughly similar in general appearance, except 
for the strong feature in the relativistic spectrum at u> ~ 0.5M -1 . Of particular interest in Fig. ||(b) is the excitation 
of the / mode (shown in detail in the inset). Perhaps the most important feature of our model is that we can compare 
the excitation given by a relativistic and a Newtonian computation. 

In order to illustrate how the methods of this paper work for a very weakly relativistic model, Fig. [| shows two 
spectra for a HIF star with R = 20M . In both, the Newtonian star-only spectrum is compared with the relativistic 
star-only (i.e. , subtracted) spectrum, and for both -R s hcii = 22M . Figure |3|(a) shows the case for a Gaussian parameter 
a = 0.01A/ -2 , while (b) shows the case a — 0.00001A/ -2 . The Newtonian and relativistic computations of / mode 
excitation (shown in the inset) agree quite well in Fig. |^(a), and there is general agreement in the overall shape of 
the spectrum but, as in Fig. 0(b), the relativistic case has structure that is missing in the Newtonian case. This 
is due to the fact that the shell stress energy source has a timescale (~ a^ 1 / 2 = 10AT) that is less than the light 
travel time across the star, and the star is not radiating in phase. By comparison, in Fig. ||(b) the shell timescale is 
~ aT 1 ! 2 ~ 300M and the slow approximation is justified. (Note that the / mode oscillations, at to ~ 0.01M -1 are 
also slow compared to the light travel time across the star.) Because of this, the Newtonian and relativistic spectra 
in Fig. 0(b) are in excellent overall agreement. 
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FIG. 1. A waveform and spectrum for an extremely compact HIF model with R — 2.5M. For the waveform in (a) the 
shell radius was taken to be large, -Rshcii = UOM, and the time scale for the shell stress energy small, a = 0.1M~ 2 , in order 
to show a clear separation of the initial shell radiation and the radition due to the stellar fluid. Waveforms are shown for the 
relativistic star+shell, and for the relativistic shell-only model discussed in the text. For the spectra in (b) the shell radius is 
small, -Rshcii = 5M, and the shell time scale is slow, a — 0.001A/ -2 , so that the model is approximately a slow motion source. 
The spectrum of the relativistic star+shell model is compared with the Newtonian shell+star model. The agreement is good, 
but is dominated by the shell radiation. For the R — 2.5M model the relativistic / mode is (0.18 + «0. 000037) A/ -1 , and the 
Newtonian / mode is (0.226 + iO. 00131) Af~ ; the least damped w mode is (0.42 + i0.022)Af _1 , while the least damped w mode 
for the shell-only model is (0.44 + i0.029) A/ -1 . 




-100 100 200 300 400 0.25 0.5 0.75 1 

u/M ci)M 

FIG. 2. Results for a HIF model with R = 5A/, i? s heii = 10M, and Gaussian parameter is a = 0.01A'/" 2 . In (a) the 
radiation waveform from the relativistic shell+star is compared to that from the relativistic shell-only to demonstrate the 
value of subtraction. In (b) the spectrum for the star-only (i.e. , subtracted) relativistic computation is compared with the 
Newtonian star-only spectrum. The / mode frequency for the R — 5M model is (0.08 + iO. 000033) Af" 1 in relativity and 
(0.08 + i0.000082)A/ _1 in the Newtonian computation. The least damped star+shell w mode frequency is (0.50 + i0.31)Af~ . 




FIG. 3. A comparison of Newtonian and relativistic star-only computations for a weakly relativistic HIF model. In both 
figures the stellar radius is R = 20M, and the shell radius is i? s heii = 22Af . In (a) the Gaussian parameter is a — 0.01M" 2 , so the 
shell is not a "slow motion" source. In (b) the Gaussian parameter is a — 0.00001M -2 , so that the shell is a slow motion source, 
and the Newtonian and relativistic results agree in all features. The / mode frequency for this model is (0.01 + i0.32)10 -6 M _1 
in Newtonian theory and approximately the same in relativistic computations. 

VI. CONCLUSIONS 

We have presented a method of probing the gravitational wave properties of neutron stars by using time varying 
stress energy in a spherical shell. In particular, we have shown that this method can give more useful answers about 
the neutron star physics than those given by studies of the scattering of gravitational waves j?J or by the response of 
the star to a close particle orbit [p|-|lO[|. The main motivation for considering such a probe is to compare relativistic 
and Newtonian computations of gravitational radiation. We have shown that the shell probe is well suited for this 
purpose, since both the Newtonian and relativistic computations can be carried out for the model. 
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Two additional features of the shell probe have been shown to be important or useful. One is the possibility of 
approximately distinguishing the radiation that can be ascribed to the star from the radiation due to the shell. Since 
the radiation from the shell is stronger than that from the the stellar fluid, this separation is valuable in bringing out 
the physical radiation that is of primary interest. 

A fundamentally important feature of the shell probe model, is the ability to choose the timescale of the shell stress 
energy. This has allowed us to direct attention to the fact that the Newtonian approximation is not only a weak 
field approximation, but a "slow" approximation. That is, the quadrupole approximation used in a quasi-Newtonian 
gravitational wave calculation supposes that the light travel time across the source is much smaller that the period 
of the waves generated. We have shown that even for a weakly relativistic stellar model there is not good agreement 
in the details of the Newtonian and relativistic computations if the timescale for the shell stress energy is short. 

The results shown in the previous section make this clear. The structure and dynamics of the R = 20M model 
of Fig. ||(a) is well described by Newtonian physics, but there are large differences between the Newtonian and 
relativistic spectra when the star is excited by a short time scale perturbation. The difference between the Newtonian 
and relativistic results is even larger for the R = 5M model of Fig. ||(b). For smaller Gaussian parameters a (i.e. , for 
driving perturbations with longer time scales) the "relativistic-only" structure seen in Figs. |(b) and |(a) deer eases. 

What becomes clear from these results is that the question of whether Newtonian physics is adequate for neutron 
star dynamics is inseparable from the question of the timescale of the excitation of the neutron star. If the timescale 
is imposed from a distance many times the neutron stars radius, then the excitation will be slow and our results 
(based on a very limited exploration of models) suggest that Newtonian physics will suffice. On the other hand, rapid 
processes, due to impacts, collapse, etc., may be have a timescale only several times GM/c 3 , the slow approximation 
may be violated, and Newtonian calculations may be significantly in error. 

As explained in Sec. I, the motivation for the shell probe for neutron star oscillations is a sequence of two questions. 
In Ref. Andersson and Kokkotas questioned whether Newtonian physics was adequate for neutron star physics. 
The second question is whether this can be adequately studied with the particle-scattering model [p-10 and its 
inflexible timescales. The results presented here suggest that in cases in which neutron stars are excited on a very 
short time scale, those particle-scattering computations are not a sufficient basis for the conclusion that relativistic 
effects arc unimportant in neutron star models of gravitational radiation sources. 

Our main purpose here has been to introduce the shell probe, and the motivation for it. We have applied it only to 
a single simple neutron star model. It is quite possible, of course, that for some equations of state the importance of 
relativistic effects might be quite different. If other equations of state are to be studied towards this end, we suggest 
that the method of time varying stress energy in a shell be considered as a good way of getting the clearest comparison 
of Newtonian and relativistc predictions. 
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APPENDIX A: THE NEWTONIAN LIMIT OF A HIF 



1. The newtonian perturbation equations and their solution 



We consider a HIF, excited by a spherical shell whose mass is much smaller than the star's mass. Assuming that the 
motion of the fluid consists of small perturbations around a spherically symmetric and static fluid ball (the equilibrium 
star), we can decompose all perturbative scalar quantities, namely the gravitational potential, the pressure, and the 
density perturbations, in scalar spherical harmonics, and can decompose the fluid's velocity in even parity vector 
spherical harmonics. We write these decompositions as 

U = U eq (r) + J2 su i(r,t)Y l0 (6) (Al) 

p = Peq (r)+ Peq Y,5h l {r,t)Y la {0) (A2) 
i 
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shcllj 



Wi(r,t) 



p eq 5[r-R]}Y w (d) 



» r = -^£|Wo(r,*)«o(«) 
i 



(A3) 
(A4) 
(A5) 



where U eq ,p eqi p eq = const are the gravitational potential, pressure and constant density of the spherical (equilibrium) 
star. The perturbation of the density has two contributions: one from the matter in the shell and the other from the 
star itself. In writing the first one, we supposed that <5Too. shell ~ <5/9shcii is dominant over all the other components of 
the stress energy tensor of the shell (weak field, slow motion approximations) and that (5Xbo jS heii is given by Eqs. (|^) 
and (|ll| ) . The density due to the fluid perturbation is zero everywhere inside the star (since the star is incompressible) , 
but not at the unperturbed surface of the star In writing the perturbation of the pressure, we used the definition 
Eq. ([It]) and the fact that p eq <§; p e q in the weak field limit. 

We then substitute these expressions in the fluid equations for a perfect barotropic fluid, which are Poisson's 
equation for the gravitational potential, 



the continuity equation, 



and Eulcr's equation, 



W l U = -Airp , 



P,t = -V.(/w) , 



(A6) 



(A7) 



v t + (v.V)v = VU Vp 

P 



(A8) 



and we keep terms only to first order in the perturbations. The resulting linearized fluid equations can then be 
reduced to two equations: one for SUi and another for Shi. These equations are the Newtonian limit of the relativistic 
even parity perturbation equations derived by Lindblom et al. pjj ]. In the special case of a HIF, they reduce to two 
decoupled second order equations, 
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(A9) 



and 
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(A10) 



Since the left hand sides of Eqs. (A9) and (A10) are equivalent to the multipole decomposition of the Laplacian, it is 
straightforward to write down the solutions that are well behaved at the center of the star and that vanish at infinity: 
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shell 



and 



Sh l (r,t)= m(ty, r<R. 



(A12) 



The three constants can be easily determined, by requiring the vanishing of the Lagrangian pressure at r = R, as in 
Eq. (|4li|), by requiring the continuity of SUi at the surface of the star, and by integrating Eq. (A9) about r — R and 
by using the linearized Eulcr equation for the radial component of the fluid's velocity. The result is 
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The frequency loi is the / mode of vibration of the star. 



(A13a) 
(A13b) 

(Al3c) 
(A13d) 



2. Damping of the / mode oscillation 

Once set into vibration at its / mode frequency by the shell perturbation, the star would oscillate forever, since 
there is no damping mechanism for the perfect fluid in Newtonian theory. This is clear from Eq. (A13c), which is 
the equation of an undamped harmonic oscillator. But in practice the star will radiate away its vibrational energy 
through gravitational waves, over a long period of time. The time, 77, for gravitational wave damping of the / mode 
oscillation can be computed, in the weak field, slow motion, approximation using energy conservation E2] p3| to be, 



41(1 - l) 2 (2l + l)[(2l - 1)\\} 2 
3(l + l)(l + 2)u 2l+2 R 2l + 1 



To introduce this damping we replace Eq. (Al3c) by 
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(For a similar procedure see Q and ^|). The retarded solution of this equation is 
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From Eq. (Al3a) and the linearized Euler equation we obtain 
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Combined with Eq. ( |A16| ) and Eq. (|A3|) for I = 2, this leads to the following expression for the quadrupole perturbation 
of the density 
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(A18) 



3. The energy radiated by a vibrating, semi-Newtonian HIF 

The energy radiated in gravitational waves by an oscillating Newtonian star can be computed by regarding the star's 
gravitational field as a small perturbation of Minkowski spacetime. For details see Ref. [Q . The energy radiated in 
the quadrupole is 



dE, 



1 



where A20 is the Fourier amplitude of the mass quadrupole, 

POO 

A 20 (oj) = / due iuu I 20 (u) . 



(A19) 



(A20) 
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In the slow motion approximation, l2o(u) is simply 



ho(u) 
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6p2(u, r)r 4 dr . 



In this approximation, Eqs. (A18) and (A20) can be combined to give the Fourier amplitude 
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(A22) 



which can then be used in Eq. (A19) to compute the Newtonian quadrupole energy spectrum. The first term in 
curly brackets in the Fourier amplitude Eq. (A22) is the contribution only of the shell to the total energy radiated. 



The second term is the contribution only of the star's oscillations to the energy. This should be contrasted with the 
relativistic procedure in which no such exact identification of the shell and star contributions can be made. 
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